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CHAPTER  I 


THE  ROOT-SQCARING  METHOD 


Historical  backgroimd 

A method  for  approximating  all  of  the  roots  of  a polynomial  equation 
simultaneously  by  using  a "root-squaring"  process  to  separate  the  roots 
has  been  variously  attributed  to  three  men;  Dandelin,  Graeffe  or  Lobachevs- 
kii  [U].  The  first,  G.  Dandelin,  in  four  appendices  to  a paper  [7]  pub- 
lished in  1826,  described  the  basic  procedure  for  "root-squaring"  and 
indicated  how  it  might  be  used.  However,  his  work  apparently  received 
little  notice  and  it  remained  for  C.  H.  Graeffe  [10]  in  I857  to  develop 
the  method  essentially  as  it  is  used  today.  In  much  of  the  subsequent 
literature  on  the  subject,  the  method  bears  the  latter's  name;  "Graeffe 
Root-Squaring."  In  I834,  lobachevskii  included  a brief  section  in  his 
book  ;^gebra  [I3]  on  the  root-squaring  technique,  and  some  authors,  par- 
ticularly Russians,  are  inclined  to  give  him  the  credit.  But,  according 
to  C.  A.  Hutchinson  [12],  the  fundamental  idea  involved  goes  back  to  two 

papers  by  Edward  Waring,  "Miscellanea  analytica,"  1762  and  "Meditationes 
analyticea, " I776. 


The  root-squaring  algorithm 

The  root-squaring  process  itself  is  based  on  the  very  elementary 
product  of  a sum  and  a difference  of  like  quantities. 


The  general  polynomial  of  degree  n with  real  coefficients, 
P(x)  = + • • . + 


n-1 


X + a 


1 


2 


written  in  factored  form  becomes 

(1.2)  P(x)  = a^(x  - x^)(x  - Xg)  • • • (x  - x^), 

with  x^,  Xg^  * • • , x^  the  n roots  of  P(x)  = 0.  Clearly,  using  (1.2), 

(’1)^  P(“3c)  P(x)  = (-1)“  a^  [(-x-x^)(-x-X2)  • • • (-x-x^)] 

[(x-x^)(x-xg)  • • • (x-x^)] 

[(x+x^)(x+xg)  • * • (x+x^)]  [(x-x^)(x-X2) 

• • • (x-x  )] 

(x^-x^)(x^-x^)  • • • (x^-X^) 

t 

(y-x^)(y-x^)  • • • (y-x^)  = Pg(y), 

and  P^(y)  = 0 is  an  eq^tion  whose  roots  are  the  squares  of  the  roots 
of  P(x)  = 0. 

The  coefficients  of 


(1-3) 


= a 


= a 


= a 


(1.4)  Pg(y)  = A y^  + A^y^"^  + * • * + A ,y  + A 
ox  n-l*'  n 

are  compactly  represented  by  the  convolution  product  [9]* 

(1-5)  (\)  = ^ (k  > n,  = 0). 

Either  from  this  or  directly  from 

(1.6)  (-if  P(-x)  P(x)  = H-  t • • • + + A^, 

the  A’s  are  foimd  to  be 

(l.T)  A,,  = (-if  iB^  + (-l)J 

*The  convolution  product /b^\  = fA  * L\  defines  a sequence 
terms  of  two  given  seguenbes  f a.  \md  /b  \}  k = 0,  1,  2,  • • • by 

\ “ %»k  - “l4-l  ^ t a,b^. 
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The  latter  form  is  the  basis  for  the  usual  organization  of  the  work  for 
a desk  computer,  n + 1 columns  are  headed  by  a^,  a^,  ‘ ‘ t Below 

each  a^  is  recorded  the  corresponding  found  by  sq'uaring  the  a^  di- 
rectly above  and  totaiUing  successively  with  alternating  signs  twice  the 
products  of  the  aj^_^  and 

until  either  a^  or  a^  is  reached.  For  a fourth  degree  equation,  for 
example: 


symmetriceOly  located  on  either  side  of  a^^^ 


^o  \-2v2  % 

k I 

Note  that  the  (-1)  in  (1.7")  is  dropped,  so  the  coefficients,  A^,  correspond 

to  an  equation  whose  roots  are  the  negatives  of  the  squares  of  the  roots 

of  P(x)  = 0,  the  so-called  "Encke"  roots.  If  this  procedure  is  repeated 

t 

on  the  A^,  etc.  m times,  the  coefficients  of  an  equation  whose  roots  are 
the  negative  2“  th  powers  of  the  roots  of  P(x)  = 0 are  obtained. 


Root  approximation 

As  suggested  previously,  the  purpose  of  root-squaring  is  to  further 
separate  the  roots  of  an  algebraic  equation.  That  is,  if  the  moduli  of 
the  roots  of  P(x)  = 0 are  unequal,  the  moduli  of  high  powers  of  these 
roots  will  differ  much  more.  The  squaring  process  is  continued  until 
the  ratio  of  the  modulus  of  any  smaller  root  of  the  derived  equation  to 
that  of  any  larger  root  is  negligible  compared  to  unity.  Suppose  the 
moduli  are  ordered, 

IxJ  > Ix^l  > • • • > )x^_J  >(x^(  , 

and  that  the  above  condition  obtains  at  the  m^  root-squaring.  Letting 
p = 2 ^ this  implies 


k 


(1.8) 


J = 1,  2,  3, 


n-1;  r > 0. 


If  Pp(x)  - 0,  the  equation  whose  roots  are  the  p—  powers  of  P(x)  = 0» 
is  divided  through  hy  a^,  the  resulting  equation  has  the  form 

(1*9)  + • • • + B tX  + B =0. 

J.  ei  n-1  n 

Expressing  these  coefficients  in  terms  of  the  elementary  symmetric  functions 

p 

of  the  roots  of  P^Cx)  = 0 and  letting  x^  = j = 1,  2,  • • • , n gives 

[14,  pp.  216-217] 


(1.10)  Bj_  = -(X3_  + Xg  + • ‘ + x^)  = -X^|l  + ^ 4-  i + . . .4-^1 


B2  = 4-  Xj_X^  + 


1 + ^ + 


= -(X^X^X  + Xj^XgXi^  4- 


= -% 


4^^ 


% * + Vn  * ==3^ 


+ + — — 4-  * 

X2  X^^ 


X 
_r 


■ V A 


X tX 
n-1  n 


+ ^ + 


% 

■ • • ♦ Wn-A) 
V2Wnl 


\ • • • \ • 

Using  (1.8),  these  are  approximated  by 
(l.li)  = -X^,  = X^)^,  Bj  = 

• • • ®n-l  = • • • X^.^,  B^  = (-1)\J^  • • • 
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Tnus , 


(1.12)  ‘ -h’i  - -h’  " ' ‘ 

J.  2 ^ n-1 


+Vi 

estimate  tlie  roots  of  P (x)  = 0 and,  hence,  the  p — powers  of  the  roots 

P 

of  P(x)  = 0. 

The  usual  criterion  for  the  temination  of  the  root-squaring  process 
is  for  certain  corresponding  coefficients  at  one  stage  to  increase,  within 
the  limits  of  precision  required,  by  squaring  at  the  next  stsige.  That  such 
"pivotal"  coefficients  exist,  regardless  of  the  number  of  roots  with  the 
same  absolute  value,  is  sho^m  by  Erwin  H.  Bareiss  [2].  He  also  shows  that 
ratios  of  consecutive  pivotal  coefficients  give  the  moduli  of  all  the 
roots . 


The  problems  of  complex  roots,  repeated  roots,  and  roots  with  equal 
moduli  are  treated  as  special  cases  in  most  general  works.  C.  A.  Hut- 
chinson [12]  summarizes  a number  of  the  most  common  types  in  terms  of 
the  behavior  of  coefficients  at  successive  stages  of  the  root-squaring 
and  the  kind  of  roots  implied.  L.  L.  Cronvich  [6]  refines  Hutchinson's 
list  in  detail  and  extends  it  to  more  cases.  Once  the  moduli  and  types 
of  roots  are  known,  the  signs  of  the  real  roots  are  easily  found  by  sub- 
stitution in  the  original  equation.  Determination  of  the  complex  roots 
is  largely  a matter  of  the  ingenuity  of  the  calculator— he  may  reduce 
the  original  equation  with  the  determined  real  roots,  utilize  the  sym- 
metric functions  of  the  roots  or  employ  various  transformations  and 
other  algebraic  devices. 
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Method  of  Brodetsky  and  Smeal 

Brodetsky  and  Smeal  [5]  avoid  many  of  the  difficulties  of  extracting 
high-order  roots  to  obtain  the  real  roots  and  of  determining  amplitudes 
of  imaginary  roots  by  applying  the  root-sq^uaring  procedure  to  P(x  + h)  as 
veil  as  to  P(x).  Assuming  h to  be  of  small  enough  magnitude  that  higher 
than  first-degree  terms  in  h are  negligible,  they  develop  formulae  for 
the  coefficients  of  P^Cx  + h)  in  terms  of  those  of  PgCx),  h and  a set  of 
^ ~ The  latter,  b^,  are,  in  turn,  defined  in  terms 

of  the  coefficients  of  P(x).  As  with  the  regular  root-squaring  process, 
this  procedure  can  be  iterated  as  many  times  as  necessary.  Finally,  then, 
the  real  roots,  with  no  ambiguity  of  sign,  and  the  amplitudes  of  complex 
roots  for  a nianber  of  cases  can  be  found  from  simple  functions  of  ratios 

of  the  b*s  and  the  coefficients  of  P (x). 

P 

Besultant  procedure 

Quite  recently — i960 — Ervin  H.  Bareiss  [2]  derived  criteria  for 
the  number  of  root-squarings  necessary  to  distinguish,  as  different, 
roots  with  moduli  in  a given  ratio.  He  also  found  corresponding  criteria 
for  the  acceptance  or  rejection  of  coefficients  in  P (x)  as  pivotal; 
i.e.,  such  that  their  ratios  approximate  the  powers  of  the  moduli 
of  the  roots  of  P(x)  = 0.  Once  the  set  of  pivotal  coefficients  and, 
hence,  the  set  of  differing  moduli,  is  determined,  the  actual  calculation 
of  roots  is  accomplished  by  a "resultant  procedure."  For  any  given 
modulus,  p,  of  any  multiplicity  for  P(x)  = 0,  the  resultant,  P(p),  of 
P(x)  and  q(x)  = x^  + px  + is  calculated.  In  general,  R(p),  is  a 
polynomial  of  degree  n in  p.  A necessary  and  sufficient  condition 
for  P(x)  and  q(x)  to  have  a common  linear  factor  is  for  their  resultant 
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to  vanish;  thus,  the  solution  of  R(p)  = 0 will  yield  all  such  factors. 

Since  complex  roots  of  P(x)  = 0 always  occur  in  conjugate  pairs,  a and 

a,  if  X - a is  a common  factor  of  P(x)  and  Q(x),  so  also  is  x - Qt.  Then 

q(x)  = (x  - a)(x  - a)  = X -(o:  + a)x  + oa  and  p = -(o:  + a)  is  real.  Also, 

2 ^ 

for  roots  of  modulus  p,  only  those  roots  of  R(p)  = 0 with  p = need 
be  considered.  Hence,  all  and  only  those  real  roots  of  R(p)  = 0 satisfy- 
ing  p = 4p  , substituted  into  q(x)  = 0,  will  yield  the  roots  of  P(x)  = 0 
with  modulus  p. 

Root-cubing 

A.  C.  Aitken  [1],  in  a brief  note,  describes  a method  of  root- 
cubing;  i.e.,  a method  for  finding  P^(x)  = 0,  the  equation  whose  roots 
are  the  cubes  of  che  roots  of  P(x)  =0.  He  shows  how,  in  two  stages, 
instead  of  the  one  stage  of  root- squaring,  the  coefficients  of  P^(x) 
can  be  determined.  The  calculations  are  similar  to  those  for  root- 
squaring in  that  they  involve  algebraic  stuns  of  binary  products  only; 
however,  more  terms  are  needed  and  the  signs  of  the  products  are  not 
simply  alternating  as  in  the  latter  process.  He  points  out  that  the 
only  significant  advantage  over  root-squaring  is  that  real  roots  are 
obtained  without  ambiguity  of  sign. 

Aitken  [l6]  has  also  extended  the  Bernoulli  method  for  the  largest 
root  of  an  algebraic  equation  to  all  the  roots.  He  employs  functions 
related  to  the  Newton  symmetric  functions  and  obtains  the  unpowered 
roots,  but  in  a somewhat  different  form  from  that  proposed  in  the 
following  chapters. 


CHAPTER  II 


GENERAL  ROOT-POWERING 

Introduction 

In  this  chapter,  four  methods  for  raising  the  roots  of  a polyno- 
mial equation  to  any  power  are  developed.  The  first  utilizes  a convo- 
lution product  involving  complex  roots  of  \inity.  The  second  is  based 
on  a special  property  of  circulant  matrices.  Ordinary  characteristic- 
number-powering  of  matrices  is  adapted  to  the  special  case  of  polyno- 
mial equations  for  a third  method,  while  in  the  fourth,  it  is  pointed 
out  how  the  relations  developed  by  Newton  between  certain  symmetric 
functions  of  the  roots  and  the  coefficients  of  P(x)  can  accomplish  the 
same  end. 

The  fundamental  theorem 

This  theorem  is  based  on  the  very  elementary  fact  that  the  p 
distinct  roots  of  the  equation,  x^-r^  = 0,  ra  constant  and  p any 
positive  integer,  can  be  expressed  in  terms  of  r and  p p^  roots  of 
unity: 

jr|  ^ ~ k = 0,  1,  2,  • • • , p-i 

2m 

where  = e ^ . The  basic  theory  of  root-squaring  and  root-cubing 
discussed  in  Chapter  I,  as  well  as  the  theory  of  the  first  two  methods 

of  general  root-powering  below,  derive  more  or  less  directly  from  the 
theorem. 
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Fundamental  theorem:  If 

(2.1)  P(x)  = a + a-x^”^  + ...+a  ,x+a  ;a,a  / O 

^ ' o ^ n-1  n ' n 

= a^(x-x^)(x-Xg)  . . . (x-x^), 

then 

(2.2)  n P(0J^X^)  = (-l)^^^'*'^^a^  (x-x^)(x-x|)  . . . (x-x^). 

j=0  ^ o X n 

That  is,  an  equation,  Pp(x)  = 0,  whose  roots  are  the  p — powers  of  the 
roots  of  P(x)  = 0 is  given  by 


n-l 


2m 


IT  P(cD'^x^)  = 0,  where  „ 

j=0  P ' P 

Proof: 


p 

m = e 


p-1  i p-1  i - 

n P(cD^x^)  = n a (m'^x^-x^  )(o)'^x^-x_)  . . 
j=0  0=0  ° ^ P 2^ 


. (oi'^X^-X  ) 
' p n' 


p-1  n - 
= IT  a n (fa^yP-x,  ) 
j=0  °k=l  P ^ 


n ^ IL 

= a^IT  IT  (oi'^x^-x,  ) 
Vi  j=o  P ^ 

= a^  n n 

° k=l  j=0 


IJ 


n p-1  /X, 


= (-1)-.^  f-l-J 

Vx^ 


k=l 


n 


= aP  IT  (-1)I>  X 
° k=l 


n 

P 

1 

- 1 

L ypi 
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= a^ 


I n (-Ip  - X)  = a»  n 

k=l  ^ ° k=l  ^ 

= (-1)“^^'*‘^V  (x-x^)(x-:!^)  . . . (x-:^). 


(2.5)  (-1) 
then 


Corollary  1:  If 

n(p+l) 


,np-l 


n P(co"s)  = As^  + as--+...+a  tS+A  , 
j=0  ^ o T.  np-1  np^ 


(2-4)  * . . . » 

That  is,  the  seq,uence  is  generated  by  the  p-fold  convolution  of 

sequences  formed  from  the  coefficients  of  P(x)  and  powers  of  co 

P 

Proof: 

P"1  . p , 

n P(co^s)  = H P (aPs)  , j = p-h. 


J=0 


h=l 


-h 


P-1 


= n P (co’^^s)  = n P (aj"^s)  , 


h=l 


h=0 


since  = 1.  For  any  particular  h, 


-h(n-l)  n-1 


1 P 


n-1  p 


= 0) 


-hn 


(%s^  + cOp  a^s^"^  + + . . . 


p n-1  p n 


-hn  / X 

- “p  ' 


la  which  pj^(s)  can  be  considered  as  a generating  function  of  the  sequence 
/ hh  I 

^ '^P-l  , P-1  ^ p-1  p-i 

II  p(OpS)  = n 0)-  Pj^(s)  = II  o-"”  n p (s) . 

J-t*  h-0  h=0  h=0 
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Now, 

P-1 

rr 

h=0 


^.^P(Pri). 

-hn  -n'^E  n 2 

n 0)  = CO  o = CO 

P P P 


If  p is  odd> 

J^P(P-I)  _ n(p-l) 

^ ^ = ^)  ^ =1  , 


_ Q(P“1) 
since  2 


is  an  integer  and  co^  = 1. 


Also,  if  p is  odd. 


P 


h=0  ■ 

p 

5 

CO 

P 

= -1,  so  that 

-hn 

jVW 

1 

tf 

»5' 

1 

H 

h=0 


Thus,  in  either  case. 


h=0 


and 


P-1 

H 

J=0 


P-1 

n 

h=0 


pi  1 

H P(o)fs)  = ri"  It  p (a)  = £ P Cs) 

h=0  *“  h=^ 


or 


(-1) 


,n(p+l)  ^ ^ ^ 


j=0  ^ h=0 

A theorem  on  convolutions  [9,  p.  251]  states  that  the  product  of  the 
generating  functions  of  two  (or  more)  sequences  gives  the  generating 
function  of  their  convolution.  Applying  this  theorem  to  the  last  equation 
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above  and  noting  that  the  right-hand  side  of  (2.3)  is  just  the  generating 
function  of^A^  , we  have  the  desired  result>  (2.4). 

Corollary  2;  = . . . = ^ = 0, 


^+1  " “^p+2  = • • • 


(2.5) 


^(n-l)p+l  “ ■^(n-l)p+2 
A^  = 0,  k > np. 


= A , = 0, 
np-1  * 


Or,  equivalently. 


(2.6)  (-l)^^P'^^)  n P(co^s)  = n p (s)  = A s^^  + A 

j=0  ^ j=0  J ° ^ 

+ . . . + A . 

np 

Proof:  From  the  fundamental  theorem,  the  right-hand  side  of  (2.2) 

will  involve  only  integral  powers  of  x.  Thus,  (2.3)  will  involve  only 
pov7ers  of  s that  are  multiples  of  p,  all  other  coefficients  vanishing. 
For  an  example,  let  n=4,  p=3,  so 


P-1 


P(x)  = a x^  + a-,x^  + 


Then, 


P^(x)  = P(x^)  P(co^x^)  P(o|x^) 

= + A^x^  + A^x^  + AgX  + 

and,  dropping  the  subscript  on  (jo. 


{^}  " {\)  * (“\)  * {-\}  ■ 


Now,  making  use  of 

85s 

03  = 03  = 03 


03  = 03  = 03,  03  = 03  =1 


and  03  + 03  = -1,  we  have 
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(“V  = 


o> 


'2^ 


ana., 


2 2 
0)  a^,  oia^,  a.y  oi  a^^. 


{‘"^\)  = %> 

I^Uins 


gives 


/ 


k 

0 

1 
2 

3 

k 

5 

6 

7 

8 


a 


a a (oi  -Ko)  = -a  a^ 
o 1'  o 1 

2 2 2 
a^a^Cdi  -Ko)  + a^  = a^  - a^a^ 

2a^a^  + a^a^  (okd^  ) = 2a^a^  - a^a^ 

+ a^a^(o)^+ai)  + a|  = a|  - a^a^  - a^a^ 

2a^a^  + a^a^Co)^ -ko)  = 2a^a^  - a^a^ 

a^aj^Cco^ -ko)  + a|  = a^  - a^a^ 

2 

a^aj^Coj  -Ku)  = -a^a^ 

2 


* 


B 


gives 


k 

0 

1 

2 

3 

k 


2 2 ^ 

“ Vo  = ° 

2 2 2 2 

V2  * = 0 

• ^Ws  * ^ 

-o4  - 


^ + ^Vl^3  ■ V2  " ^ 


a 

o 

^2 

"5 


0V3 

2 2 
+ a^B.^  - a^a^  = 0 

2 2 


" ®'o®'2®3  ■ %W  V2  " + V3 

^ - a^a^a^  + 2a^a^a^  - a^a^  = 0 

5a^a^  - ^a^a^a^  + 5a^a^^  - ^a^a^a^  + a^ 


6 


14 


7 ■ V5  • VA 

2 


8 


a^a^  + 2a^a^a^  - a^a^a^  = 0 


%\  - ^ ^2^1  ’ 4%  ^ “ ^2^5  ■"  4%  - %% 


a,a,.  = 0 


9 5a^_aJ  - ^a^a^a^  + ^ 

1^4  2 2 2 2 ^ 
10  a^aj^  - a^a^^  + a^a^  - a^a^  = 0 

n 2 2 ^ 

H = 0 

12  a? 


a^a^a^ 


Circulants 


th 


The  equation  whose  roots  are  the  p — powers  of  the  roots  of  P(x)  = 0 

can  be  obtained  also  by  means  of  circulants.  Form  the  p - square 

(p  > n)  circulant  A (x)  = 

P 


(2.7; 


0 


n 


a x-“ 
o 


0 


n-1  n 

, a 
o 


a^x  ^ a,x^ 


n-1  n 

X ^ a x^  0 
o 

n-2  n-1  n 

^n-1^  • • • a^x  ^ a^x  ^ a^x^ 


hn  ^n-2^  • . . a^; 


n-3  n-2  n-1 

• • a,x  ^ a_x  ^ a,x'^ 

P 2 1 


0 


0 


0 


n 


• a x-^  0 

o 


n-1  n 
,x  ^ a x^ 


• • "1^ 

n-2 

• • a^x  ^ a^x  ^ 

n-3  n-2 

• • P a^x  ^ I 
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• • • 

• • • 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• • • 

5 It  5 

• 

• 

• 

• 

1 

• 

2 

a ^x^  a j x^  a 
n-5  n-4  n-J 

• • 0 

0 

0 

...  a -X^ 
■ n-1 

25^ 

' 1 

a ^x^  a ,x^  a , x^ 
n-2  n-3  n-4 

* • 0 

0 

0 

• • • & 

n 

a x^ 
n-1 

i 2 5 

n 

a tX^  a ^x^  a x^ 
n-1  n-2  n-5 

* * a x^ 
0 

0 

0 

...  0 

a 

n 

(2.7)  Then,  the  determinant  of  A has  the  value  [J+,  pp.  76] 


P “ ~ 2 

I A I = H [a  + a + a + 

1 p1  t n n-1  t n-2  t 


+ a^x  oj 


n-1 
p n-1 

n 


+ a x^c3^ 
o t 


*n 

1' 

n 

Jr 

= n 

t=i 

a 

0 

+ a^ 

OD.x^j  + . . . + a , 

t / n-1 

0)  x^ 
t 

+ a 

n 

}J  — 

= n P(cD  x^)  , where  o), , t = 1,  2,  . . . , p are  the  p 
t=l  ^ ^ 

distinct  p^  roots  of  unity.  Thus,  |a  (x)j  =0  represents  the  desired 
til 

a — degree  equation  by  the  fundamental  'theorem. 

Note  that  A consists  essentially  of  two  parts:  an  upper  diagonal 

Jr 

part  consisting  of  n+1  lines  of  identical  elements  parallel  to  the  main 

diagonal  of  a^’s,  and  an  n- square  lower  triangular  part.  We  form  the  new 

matrix,  B (x),  from  A (x)  by  multiplying  the  rows  of  the  latter  succes- 
r-i  P 

sively  by  x“?” , where  r is  the  number  of  the  row,  and  factoring  the  same 

‘^ll 

quantity  from  the  r — colinnn.  This  removes  x completely  from  the  upper 


diagonal  part  and  reduces  each  of  its  exponents  in  the  lower  triangular 
part  to  unity.  That  is. 
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(2.8) 


B (x) 

P' 


a 

S’  •> 

S 

• 

• 

• 

a. 

a 

0 

• 

n 

n-1 

n-2 

1 

0 

0 

a 

n 

s 1 
n-1 

• 

• 

• 

ag 

a 

0 

• 

• 

• 

0 

0 

a 

n 

• 

• 

• 

3-2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

0 

0 

• 

• 

• 

• 

• 

• 

a 

0 

a X 
0 

0 

• 

• 

• 

• 

• 

^2 

a^x 

a X 
o 

0 

• 

• 

• 

• 

• 

• 

^2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a ,x 
n-3 

a 1 X 
n-4 

a 

n-5 

• 

• 

• 

• 

• 

• 

s , 
n-1 

^n-2 

a „x 
n-2 

a ,x 
n-3 

a I X 
n-4 

• 

• 

• 

0 

« • 

• 

• 

a 

n 

\-l 

a -X 
n-1 

a ryX. 
n-2 

n-3 

• 

• 

• 

a X 
0 

0 

• 

• 

• 

0 

For  example,  with  n = 3,  p = 8, 


A3(x)  ^ 


i 2 3 

^ H H 

a^x  a^x  a^x  0 

12  3 
H . . H H 


0 


a^x  a^x  a^x  0 


0 


0 

0 


5 

a X 0 0 0 

o 

2 3 

H H 

a^x  a^x  0 0 

12  3 

H H H „ 

V V V ° 


0 

12  3 

H H H „ 

a^x  a^x  ax  0 

1 o 

i 2 3 

H 8 H ^ 

a^x  a^x  a^x  0 

12  3 

H H 8 

a^x  a_x 


2 

0 


a^x 


0 


o 

1 2 

8 8 

a^x  a^x 

1 

8 

^3  V 
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and 


BgCx)  = 


h 

^2 

a 

o 

0 

0 

0 

o\ 

0 

ag 

\ 

a 

0 

0 

0 

0 

0 

0 

a-2 

a 

0 

0 

0 

0 

0 

0 

^2 

a 

o 

0 

0 

0 

0 

0 

"•2 

a 

0 

a X 
0 

0 

0 

0 

0 

^2 

a^x 

a X 
o 

0 

0 

0 

0 

^2 

ia^x 

a^x 

a X 
0 

0 

0 

0 

0 

^5. 

Clearly,  because  of  the  nature  of  the  transformations  on  A J A 

p7  P 

=j  Bpjj  andjB^(x)j  = 0 also,  represents  an  equation  whose  roots  are  the 
til 

p — powers  of  the  roots  of  P(x)  = 0. 

By  interchanging  appropriate  columns  of  the  transpose  of  B^(x), 
we  obtain  the  new  matrix. 


(2-9) 


n 


a^x  a^x  . 

. . a -X  a 0 
n-1  n 

• 

• • 

• 

• • 

a^x  a^x  . 

n-2  n-1  n 

0 . 

• • 

• 

• • 

• • 

• • • 

• 

• 

• 

• 

• 

• • 

• • • 

• 

• 

• 

• 

• 

• • 

• • • 

• 

• 

• 

• 

• 

• 

• • 

• 

• • £L 

O 

^2  • 

. . a 

n 

0 

• • 

• 

. . 0 

a 

a,  • 

• • Ql 

0 

1 

n' 

a 0 
n 

• • 

• 

• • 

0 

a . 

o 

n- 

• • 

• * • 

• 

• 

• 

• 

• • 

• • • 

• 

• 

• 

• 

• 

• • 

• • • 

• 

• 

• 

• 

• 

. . 0 

• 

• • 

• 

• • \ 

ag  a^  . 

. . a 0 

n 

• 

• • 

• 

• * %j 
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a form  in  which  the  coefficients  of  P(x)  occur  in  their  natural  order  in 
rows  and  x occurs  only  in  the  upper  triangular  part  of  the  n-square  upper 
left  section. 

The  preceding  example  -n=3,  p=8-  becomes 

/a^x  a^x  agX  a^  0 0 0 o\ 


Cg(x)  = 


0 

a X 
0 

a^x 

^2 

^5 

0 

0 

0 

0 

0 

a X 
0 

^2 

0 

0 

0 

0 

0 

a 

0 

s-2 

0 

0 

0 

0 

0 

a 

o 

\ 

e-2 

0 

0 

0 

0 

a 

0 

\ 

&2 

&2 

0 

0 

0 

0 

a 

0 

e-2 

0 

0 

0 

0 

V 

Clearly,  here 


(2.10) 


= |cj,(x)| 


and  this  is  the  form  to  be  used  in  the  sequel.  We  employ  a Laplace  ex- 
pansion of  Cp  j on  the  first  n rows  and  consider  only  the  cases  of 
P = 2n;  for  n<p<2n,  l^.^!  would,  in  general,  be  of  low  enough 
order  that  any  convenient  method  of  expansion  could  be  used. 

With  this  condition,  p ^ 2n,  is  partitioned 


(2.11) 


pi  pt 
11  ^12 

pi  pt 
\ 21  ^221 


where  is  n by  2n.  If  p = 2n,  is  just  the  first  n rows  of  C^, 
and  C£2  and  omitted.  For  p > 2n,  will  have  orJ^  zero  ele- 

ments. Thus,  in  a laplace  expansion  on  the  first  n rows  of  C , the  set 

P 

of  non-zero  minors  on  these  n rows  is  invariant  for  all  p ^ 2n  and  these 


19 


minors  are|^j  in  number.  Hence,  the  expansion  consists  of  the  sum  of 
I jjj  binajry  products  of  an  n-square  and  a (p-n)-square  determinant  in 
which  only  the  algebraic  complements  of  the  minors  of  vary  with  dif- 
ferent p. 

For  p = 2n,  the  product  of  a minor  from  and  its  algebraic 
complement  from  and  can  be  characterized  by  a sequence  of  the 
integers  1 through  2n;  the  first  n of  them  indicating  the  columns  of 
C'  chosen  for  its  minor,  the  last  n the-  columns  of  C' , in  the  algebraic 
complementj  Cgg  constitutes  the  last  p - 2n  columns  of  every  algebraic 
complement.  The  sign  of  the  product  is  determined  from  the  number  of 
inversions  in  this  sequence,  since  this  will  correspond  in  class  to  the 
number  of  transpositions  of  columns  necessary  to  put  the  minor  of  in 
leading  position  of  C^.  Note  that  permutations  within  the  first  or  last 
n terms  of  any  given  sequence  yield  the  same  term  of  the  expansion  of 
. To  illustrate  for  n = 5,  p = 8: 

= I = 


O’  O’ 
^11  ^12 

O’  O’ 
^21  ^ 22 


a X 
' o 

0 

0 

0 

0 

a-2 


a^x 

a X 
o 

0 

0 

0 

0 


8-2^ 

0 

0 

\ 

1° 

0, 

a^x 

8-2 

0 

0 

a X 
o 

^2 

^0 

0^ 

0 

a 

o 

^2| 

f^3 

0 

0 

0 

a 

o 

^2 

^3 

0 

0 

0 

8- 

O 

^1 

8-2 

0 

0 

0 

0 

a 

0 

0 

0 

0 1 

Vo 

a 

0 
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and,  for  example. 


2^  3 ~ ” 


a^x 


V ^2 


B-r- 


0 

0 

^2 


0 

0 

0 

0 


2 

a 

o 

0 

0 


0-2 

a 


^2 

"l 


negative,  since  2,  if,  5^  1^  3^  shows  five,  an  odd  number,  of  inver- 
sions. Similarly, 


3f  2,  4,  5I0  “ + 


a X a^x  0 
o 2 


a^x  0 


ax  a^ 
o 3 


0 

0 

0 


0 

0 

0 

0 


0 

0 

0 


‘‘2 


and 


2,  3,  If;  1,  5,  6|g  = - 


a^x  a^x  a^ 

V ^1^  ^2 


Jo 


V \ 


0 

0 

0 


2 

a 

o 

0 

0 


‘"2 

a 


Now,  consider  the  matrices  for  any  two  consecutive  values  of  p, 
say,  r ^ 2n  and  r + 1.  As  established  above,  C*  is  invariant  in  C and 

JJL  p 

^r+l*  ^21  ^eros  at  the  top  of  the  first  n columns  and  at 

the  bottom  of  the  last  n columns  in  going  from  to  C while  Cpp  is 
augmented  by  a row  of  zeros  at  the  bottom  and  a new  column,  beginning 
with  a^  at  the  bottom.  The  minors  from  C^)  occurring  in  the 
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expansion  of  can,  in  turn,  be  expanded  on  elements  of  their  first 

or  last  rows;  thus,  we  compare  in  and  as'foUows; 

I.  For  expansions  on  first-row  elements. 

1.  The  first  row  of  in  has  all  zeros  except  in 

columns  n+1,  n+2,  . . . , 2n  and  2n+l.  The  subscript  on  the 
a in  each  of  these  non-zero  positions  is  n+1  less  than  the 
niomber  of  the  coliimiv. 

2.  If  the  first  row  of  in  is  deleted,  the  remain- 

ing  matrix  has  all  zeros  in  the  (n+l) — coliimn;  if  this  column 

zeros  is  removed  and  the  following  columns  each  moved  one 

to  the  left,  the  resulting  matrix  is  identical  to 

in  C . 
r 

II.  For  expansion  on  last-row  elements. 

1.  The  last  row  of  in  consists  of  zeros  except 

for  the  last  element  which  is  a^,  and  the  first  n elements 
where  the  subscript  on  the  a corresponds  to  the  number  of 
the  column. 

2.  If  the  last  row  of  is  deleted,  the  remain- 

ing matrix  has  all  zeros  in  the  n—  column;  if  this  column 
of  zeros  is  removed,  the  first  n-1  columns  moved  one  to  the 
right  and  the  last  column:  placed  in  number  one  position, 

^ the  resulting  matrix  is  identical  to  (Cl  C ) in  C . 

21  22  * r 

The  above  observations,  I and  II,  can  be  used  to  derive  the  foim  of 
recurrence  rules  for  the  expansion  of  | Cp(x)|  for  p > 2n.  Consider  any 
particular  sequence  designating  one  of  the  binary  determinant  products 
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in  C 


r+1 

(2.12) 


say 


r+1 


• • • 

where  5 n,  n < M^  = 2n  and  j + k = n.  Then,  representing  the  minors 
from  ^22^  Jnst  the  last  n terms  of  our  sequences  and  assuming 
the  N*s  and  M's  in  increasing  order,  we  expand  j . . . 

Mj^l  first  row,  using  I. 

If  Mj_  = n+1, 

;N^,  . . . ^Nj,n+1  1^,  • • • # 


r+1 


(2.13) 


(2.14) 


(-l)'^a^  1 . . . , M|^-l,  2n[^. 

If  Mj_  n+1, 

I . . . , 


= (-l)'^[aJ;Nj_,  . . . , N^, 

-aJ;Nj_,  . . . , Nj,  M^-1,J^-1, 


^Mj^-l7  2n 

, Mj^-l,2n 


+ • • • 


vk-1 


+ (-l)  " %\i\,  • • • , Kj,MJ_-l,^^-l,M^-l,  . . ., 

2n|  ^ 

+(-l)\|  ;Nj^,  • • • < Hj,  M^-1,1^-1,  . . . ,M|^-l| 

Where  s = M^-n-1,  t = J^-n-1,  . . . , z = M^^-n-1. 

This  reduction  by  elements  of  the  first  row  is  summarized  by  the  fol- 
lowing rulest 
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1.  If  n+1  appears  in  the  half -sequence^  the  reduction 
consists  of  a single  term  based  on  = n+1;  other- 
wise, each  integer,  greater  than  n yields  a term. 

2.  The  sign  of  a term  is  + or  - according  as  the  is  in 
odd  or  even  position. 

3«  The  subscript  for  a is  M^-(n+l). 

4.  The  determinant  half-sequence  for  a term  based  on 
is  obtained  by  repeating  the  integers  less  than  n+1, 
deleting  the  M^,  reducing  each  other  integer  greater 
than  n by  one  and  appending  2n  at  the  end. 

5*  If  n+1  does  not  appesir,  an  additional  term  with  factor 


(-1)  a^  and  determinant  half -sequence  having  integers 
greater  than  n reduced  by  one  is  included. 

These  rules  are  applied  to  the  preceding  examples  (n=5): 


i ; 1,  3,  6 

II 

CO 

1 ; 2,  k,  5 

03 

II 

P 

0 

1 J 5^  6 

03 

II 

i 3^  6 rj  “ 8^  I j 1,  3>  5 j Y 


) ^2  1'’  ^17  i 5I 

. . , Nj,  . . . , 


7 

r+1 


Rules  for  the  expansion  of  j ; . . 

by  elements  of  the  last  row  are  deduced  similarly  from  observations  II. 
First,  we  think  of  moving  the  last  column  of  the  determinant  into  leading 
position  and  the  last  row  into  top  position,  leaving  the  order  of  the 
remaining  rows  and  columns  unchanged.  Since  this  would  involve  an  even 
number  of  transpositions  in  total,  the  resulting  determinant  would  be 
equivalent  to  | ; . . . , K^,  . . . , | 


• > 


2k 


If  = n. 


(2.15) 


(2.16) 


I ; . . . , * ' * •»  \ I r + 1 

~ (”l)  j >1^^2^*^1^N2+X^  ...  t , , f j j. 

If  Nj  / n,  expanding  on  the  new  top  row  gives 

; N^,N2,  . . . , * • • ^ ^1  r+1 

= % I ; %+l.  Ng+l,  . . . , Nj+1,  ^ I j. 

■"  S-g  i 1^  ^2+1^  • . . f ‘ • f ^ I ^ 

+ a^  1 ; 1,  N^+1,  N^+1,  . . . , ^^  * * * ^ \ | j. 


+ (~l)  ^"2  I J 1^  N^+1^  ‘ • • > ’’  * f J, 

where  s = t = N , . . . , z = N..  This  reduction  hy  elements  of  the 
last  row  is  summarized  by  the  following  rules: 


1*  If  II  appears  in  the  half -sequence_,  the  reduction  consists 
of  a single  term  based  on  = n;  otherwise,  each  integer, 
N^,  less  than  n yields  a term. 

2.  The  sign  of  a term  is  + or  - according  as  the  is  in  even 
or  odd  position.  . 

3*  The  subscript  for  a is  N^. 

4.  The  determinant  half-sequence  for  the  term  based  on  H is 

i 

obtained  by  repeating  the  integers  greater  than  n,  deleting 
the  N^,  increasing  each  other  integer  less  than  n by  one 
and  appending  1 at  the  beginning. 
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5*  If  n does  not  appear,  an  additional  term  with  factor 
a-Q  determinant  half -sequence  having  integers  less 
than  n increased  by  one  is  included. 

These  rules  are  applied  to  the  preceding  examples  (n  = 5); 

I ; 3,  6 I g = + ^5  ; 1,  2,  6 


J 2,  4,  5 
; 1,  5,  6 


& - % 


8 


= a 


; 3,  5 

y ” ^^2  1 > 5 1 

; 2,  5,  6 

Y - 1 ; 5,  6 

7 

7 ■* 


The  most  economical  procedxire  would  be  a combination  of  the  two 
sets  of  rules: 


1.  If  n or  n+1  appears  in  the  half-sequence,  use  the  rule 
giving  only  a single  term. 

2.  If  half  or  fewer  of  the  integers  in  the  half -sequence  are 
l®ss  than  n,  expand  on  the  bottom  row  (last  set  of  rules); 
otherwise,  expand  on  the  top  row. 

If  the  determinants  of  the  minors  from  C^)  for  p = 2n  are  evaluated 
and  coded  by  their  corresponding  half-sequences,  each  higher  order  deter- 
minant can  be  evaluated  successively  by  the  above  rules  in  terms  of  the 
next  lower  order.  It  is  observed  that  the  number  of  integers  greater 
than  n in  any  half-sequence  corresponds  to  the  power  of  x associated 
with  it;  also,  that  the  same  power  of  x is  associated  with  every  half- 
sequence in  its  reduction  by  the  above  rules.  Thus,  the  coefficients 
of  Pp(x)  can  be  calculated  independently  of  each  other. 

Matrix- powering 

Without  loss  in  generality  for  our  purposes,  we  may  assume  the  coef- 
ficient, a^,  of  the  highest  power  of  x in  P(x)  to  be  one.  Thus, 


26 


The  companion  matrix  of  (2.17)  can  be  written 


(2.18) 


A 


-a 


1 Oj 


0 


0 


n-1 


-a 


0 


0 


n 


Then 


(2.19)  I A - Al  1 = (-1)“  P(A) 


and  the  characteristic  numbers  of  A are  just  the  roots  of  P(x)  = 0. 

It  is  also  known  that  the  characteristic  numbers  of  any  positive 
integral  power  of  A are  the  same  power  of  the  characteristic  numbers 
of  A.  Hence,  powering  the  companion  matrix.  A,  can  provide  a method 
of  powering  the  roots  of  P(x)  = 0. 

Clearly,  successive  squarings  of  A,  A^,  A^,  etc.  would  be  analogous 
to  the  polynomial  root-squaring  discussed  in  Chapter  I.  However,  to  be 
useful,  some  criterion  for  the  tennination  of  the  matrix- squaring  process 
and  a subsequent  method  for  calculation  of  the  separated  characteristic 
numbers  would  need  to  be  derived;  this  writer  has  been  unable  to  do  so. 

Powering  of  A might  be  used  to  obtain  equations  whose  roots  are 
consecutive  powers  of  the  roots  of  P(x)  = 0 and,  thence,  simplify  the 
calculation  of  the  original  roots  as  discussed  in  Chapter  III.  That 
is,  the  usual  polynomial  root-squaring  coiild  be  carried  out  to  the  power, 
say,  p = 2 , necessary  to  obtain  the  desired  precision.  Then,  the  matrix. 
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A,  could  be  powered  by  any  convenient  method  to  obtain  etc. 

These  consecutive  powers  of  A are  particularly  easy  to  calculate  because 


(2.20)  A^  = 


of  A. 

"u 

Suppose 

bi2  . . . 

• • • 

^2n 

• 

• 

• 

• 

• 

• 

• 

• 

• 

fnl 

o • • • 

n2 

b 

mu 

^illtiplication  of  A^  on  the  left  by  A,  (2.l8),  shows  that  each  of  the 
rows  of  A^  is  shifted  down  one  and  the  last  row  is  deleted.  By  multi- 
plication of  A^  on  the  right  by  (2.18)  the  new  first  row  is  easily  ob- 
tained as  indicated: 


/ 

Pi^n'^^12" 

a^bii+b^^  . . . 

“^j^U'*‘\,j+l  * • * 

(2.21) 

"ii 

^12  • • • 

\j  • * * 

^l,n-l 

'^In 

II 

^21 

i: 

^22  • ■ • • 

• 

• 

• 

^2j  ' • • 

• 

• 

• 

^2,n-l 

• 

• 

• 

• 

• 

• 

\ Vi,i 

^n-1,2 

^n-l,j 

^ -I 

n-l,n-l 

Vi 

Finally,  the  characteristic  equations  for  etc.  are  the  desired 

equations. 


Newton  symmetric  functions 

A final  method  to  obtain  equations  whose  roots  are  consecutive 
powers  of  the  roots  of  P(x)  = 0 that  might  be  used  in  conjunction  with 
polynomial  root-squaring  can  be  based  on  the  relations  between  the  Newton 
symmetric  functions  and  the  coefficients  of  P(x)  in  the  form  (2.17). 
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These  relations  can  be  expressed  either  in  determinant  form  [5,  p.  86] 
or  in  terms  of  a sequence  of  equations  [15,  p.  26l]  from  which  the 
determinants  are  easily  derived. 


Let 


(2.22) 


where  the  x 
Then, 


n 


are  the  roots  of  P(x) 


or,  in  general. 


2ag 


1 


0 

1 , 


(2.23)  = (-1)^ 


1 


0 

0 

0 


^ Vl  \-3  * * * ^*1 

for  k = 1,  2,  3^  ’ ' • ) ond  a.  = 0 for  j > n. 

J 

The  coefficients,  a^,  can,  in  turn,  be  expressed  in  terms  of  the  S : 


or,  in  general. 
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(2.2^) 


1 

s. 


0 

1 


0 

0 


0 

0 

0 


\ ®k-l  \-2  \-3 


for  k-1,  2,  3^=0  for  k > n.  The  sequence  of  equations 

fron  which  these  determinants  are  derived  is 
(2.25)  + a^  = 0 

^2  W + 2a^  = 0 

S ^2^1  + = 0 


/ + agSj^_2  + • • • + ^ ° > 

where  a^  = 0 for  j > n. 

+ v» 

Now,  suppose  the  equation  whose  roots  are  the  p—  powers  of  the 
roots  of  P(x)  = 0 is  written  in  the  form 

(2.26)  P^(x)  = x^  + A^x^"^  + + . . . + = 0 . 

That  is,  the  roots  of  (2.26)  are  x?,  xL  . . . , . 

•L  c n 

Then,  clearly. 


(2.2T)  = 


50 


- (-if 


n! 


0 

0 

0 


S S,  , V 

np  (n-l)p 


or,  alternatively, 

(2.28)  + \ = 0 

®2p  * Vp  + ^ = 0 

%p  * Vap  + + 3A3  = 0 

0 

®np  + V(n-l)p  - * ■ ■ ■ * Vl^p  + = 0- 

^p"  ^2p^  * * * -^  ^np  obtained  from  (2.25)  or  (2.25)  and  A^,  A^, 

/.  . . , A^  from  (2.27)  or  (2.28).  Obviously,  the  coefficients  of  P (x), 

p+l'‘ 

Pp+2^^)^  etc.,  can  be  found  similarly. 


CHAPTER  III 


ROOT  DETERMIMTION 


The  identification  of  the  types  of  roots— real,  imaginary,  equal 
moduli,  etc.— as  well  as  their  numerical  determination,  is  considerably 
simplified  by  having  a sequence  of  equations  whose  roots  are  consecutive 

powers  of  the  roots  of  P(x)  = 0.  For  the  equation  whose  roots  are  the 

■th  „ . 

powers  of  the  roots  of  P(x)  = 0,  we  will  use  the  notation 

(5.1)  P (X)  = aW  x . 0 

“ « J-  n-x  n 

Case  I.  Roots  all  real  with  unequal,  absolute  valnpg 

As  with  root-squaring,  a power,  p,  of  the  roots  can  be  found  so 
that  the  ratio  of  the  p^  power  of  any  root  to  the  p^  power  of  the  next 
numerically  larger  root  is  negligible  compared  to  one  in  absolute  value; 


(3.2) 
where 

(3.3) 


^+1 


< < 1,  k = 1,  2,  3, 


, n-1 


H > hh  ■ • • N > • • • • 


Again  using  the  Newton  elementary  symmetric  functions  on  P (x)  gives 

a(p) 

(3.4)  1 


I ^ ==!/ 


= “X, 


,(p) 


= X^X^ 


1 + -^  + • 
x^ 


x^  x^ 

• + -1 . 
^ X? 


xP  ,x^^' 

n-1  n 
x^x^ 
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32 


= 


„(P) 

— = (-l)^[xjx|  . . . + xjx|  . . . + 


+ X^  X^  ....  X^3 

n-k  n-k+1  n 


^2 


= • • • ^ 

= ...  7^ 

(-l)'^^lx^xP  . . . 

• • ■ =^+1 


4 , 

1 + + . . . + 


4 


;+l 


=4. 


:+l 


x^  , . . . x^  ' 
n-k n 

x^  . . . x^ 

1 k y 


1 , ^+2  ^ ^ n-k-1  “ * • 

•Im  T T • • # *r 


^ - ■ -^*1 1 


afe> 

n 

aP 

o 

Also, 


= (-l)°x?x^  . . . x^  . 
J-  c n 


(5.5)  n^P+1) 

^ . . . xP 


+1 


. +X^'*'^X^^^  vP+1 

n-k  n-k+1  * * * ^ 


n 


• + 


vP+1  ^+1 

^n-k  • • • ^ 

vP+1  ^1 

x^  . . . x^ 


for  k-1,  2,  . . . ,n.  Cleaurly,  then. 


JP+1)  >) 

/ T zT  \ 'k+X  • 

.(p+l).(p) 

K K+1 
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k = 0,  1,  2,  • • • , n-1; 

and  the  roots  are  found  with  no  ambiguity  of  sign. 

From  the  relations  (3*^)^  (5 *5),  sjid.  an  analogous  one  for  p + 2, 
we  note 

JP+1)  q(p+2) 

(3*7)  —7 — V = a * • • X,  = —7 ipr  • 

^ ' „(p)  o 1 2 k (p+1) 

Thus,  the  power,  p,  for  (3*2 ) to  hold  is  reached  when  the  corresponding 
coefficients  in  three  consecutive  root-powerings  are,  within  the  limits 
of  precision  required,  in  geometric  progression. 


Case  II.  Real  roots  of  equal  absolute  value 

Since  mxiltiple  roots  may  be  determined  by  other  means  (g.c.d.  of 
P(x)  and  P'(x)  ),  we  consider  only  the  case 

(5-®)  1 ==ic-i|  ^ = -Vi K4ai  • 

Again,  p can  be  found  such  that 


{3.9) 


— — r < < 1 and  ’ — — < < 1. 


J=^- 


The  symmetric  functions  involved  are 


(3.10) 


Vi 


= (-l)^"^(xP  • • • + xj  • • • + • • • + xP_. 


k+1 


. . . .,P 


^k-1 


1 + ^ 


^-1 


x^) 


x^  • • • x^' 
n-k+1 Jn 

x^  • • • 

1 ^-1 


+ • • • + 
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o(p) 

\ 


= (-1)'^  + xj  • • • 


k+2 


' + V 1 

n-k-1 


• • :^1 


= (-1)“2^  • • • =^-l  4 


V ^ 


1 + -^  + . . . + n-k-1 


^,■■■4  1 


«(p) 

\ 


= (-1)^  2:^  • • • f if  p is  even. 


= (-1)’=  • • • =^.1 4 


k+2 


1 + • • • + 


x^  • • • x^  \ 
n-k-1  n 


\ ^-1^+2/ 


= (-l)^x^  • • • P 


4^1 


- (-1^^  (xj  ■ • • xP.,xgxP^^  * x^ 


= (-l)^*^=i 


+ • • 

• + x^  , , • • 

n-k-1 

1 ^ ^+2 
1 + + • • 

k 

X? 

1 k+1 

x^  • • 

a.  n-k-1 

n"' 


• x^^ 
n 


X^^  . . . X^  1/ 

1 k-1  k k+1^ 


^ • . . xj.,xf 


■with  similar  expressions  for  p + 1 and  p + 2.  Since  successive  powers 

alternate  in  parity,  will  alternate  between  (-l)^2a“x“  * * * 

/ ,\kmm  mm 

"*  * * * ^-1^+2  ^ “ P#  P + P + 2,  etc.  Because  of 

(3.9),  the  former  for  even  order- powering  will  tend  to  increase  by  nu- 
merically larger  factors  than  the  latter  odd  order-powering.  That  is. 
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this  case  will  he  recognized  by  the  irregular  behavior  of  a^™^  in- 
creasing by  numerically  larger  factors  for  even  order-powering  than 
for  odd  order-  i)owering. 

Immediately  from  (3»10), 


,(P+1) 


,(P) 


\-i  \-i 


so  that 


Jp+1)  _(p) 

,o^  . 2 ^ ^+1  ^-1 

(5-^)  ^-TiSIT-TiT 
\-i  Vi 


Case  III.  Complex  conjugate  roots, 
Suppose  ' 

(3.15)  " Pj,  and 

Also, that  p is  such  that 


.P-1 


{5.Ht) 


K:i 


< < 1 and 


-1 

k+2 


X? 


.P-1 


< < 1 . 


The  pertinent  symmetric  functions  become 

(3.15) 


= (-1)  ■ + xj  . . . xP_2xg 


+ x^ 

n-k+1  * 


n' 


= (-i)’^-'-=i  ■■  - 4.1 


' x^ 


1 + 


4-d 


• + 


x^  , ^ 
n-k+1 


. . . x^' 

n 


x^  . . . x^  . 
1 k-1 


1 (-1)'^-^ 


, since 


=^-1 


K-d 


< < 1,  etc. 
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“ • • • =^-1=^+2 

+ • • • H"  X^  ^ • • • X I 

n-k  . nJ 


(-1)^  :^  . . . 2 cos  p ^ [l  + 


.+2 


+ • • • 


pj  2 cos  p ^ 


, . . . 

n-k  n 


X?  . . . xf  -pf  2 cos  p A 
1 k-l'^k  ^ ^ 


= (-1)^2^  . . . ^ ^ ^ > since  ^ < < 1,  etc. 


\+l 


(-1)  (xj  . . . 3^_3_P^+  . . . + • • . xJ) 


^ XJ  X^  . . . X^_Q_Pjj  • 


From  these  relations, 

(3*l6)  = -2p^  cos  p ^ , 

Vl 

a,^P) 

^+1  . 2p 

JFJ  = Pk  • 

\-l 


Similarly, 

(3.17) 


«(P“1)  (p-1) 

^ . -apP-i  C03  (p-i)  p2(P-i)  , 


®k-l 


_^(P-1)  “ ^k 
®k-l 


(5.18)  4^+1)  , 4^+1)  , ^ 

jSHT  “ (p+i>  ^ ' ^rj  = 


k-1 


^-1 
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Thus,  immediately j 


(5.19) 


(p)  (p-1)  (p+1) 

^+1  ^-1  . ^+1 
"Tp7  * "7(Ffr7 

s-k-l  ^+1  ®k-l 


a,(p) 

\-l  . 2 

^+1 


From  (3.18)  and.  (3.19) 

13^*1 

2p^  cos  (p^  + 


= -2p^“^  cos  (p^  + ^)  = 


„(P+1) 

. \ 


(3.20) 


„(P+1)  g,(p) 

\-l  \+l 

o(p+i)  ■ 

^+1  ^-1 

(p+i)J:p)  * 

^+1  \-i 


Now, using  the  relation  cos  (p^i  -<f>)  + cos  (p^  + ^ ) = 2 cos  p^  cos  ^ 
together  with  (3.17)  and  (3. 20), we  find  that 

(3.21)  -4p^-^  cos  cos  # i . 

®k-l  \+l  Vl 

f ♦ 

Finally,  using  (3.I6),  (3.19)  and  (3. 21), we  obtain 

^P-l)  a(P^^)3(^ 

HFi)  TTp+rrtp7 

^ ^+1  ^-1 

cos  ^ 


^ =i!i 


(pm 


1 

rjp+Djp)-] 

\+i  Vl 

a 

_(p+i)jp) 
LVi  ViJ 

Consideration  of  (3.15)  and  the  corresponding  expressions  for  p-1 
and  p+1  indicate  that  the  coefficients  a^“^  and  a^“|  will  show  behavior 
Uke  case  I for  m = large  pj  i.e.,  a^^^,  and 

®k+l^^  will  become  geometric  progressions.  However,  will 
show  irregular  behavior,  both  in  sign  and  magnitude,  due  to  the  factor 
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cos  TMs  irregularitjj  coupled  with  the  regularities  of  the  adjacent 
coefficients,  provides  the  criterion  for  the  recognition  of  the  pair  of 
complex  conjugate  roots. 

Case  TV.  Roots  with  the  same  modulus 

Suppose  several,  U,  roots  have  the  same  modulus,  hut  different 


amplitudes,  That  is, 

i^ 


'^k+r 


, r - 0,  1,  2,  . . . , 


(5-25)  Vr  = 

and 

(3.21*) 

The  dominant  parts  of  the  pertinent  symmetric  functions  for  large  p are, 
then, 

a^(p) 

(3.25)  ^ 4 


4''^  v-i 

— S cos  p 


4^’ 

x+1  . 


= Z 3_ 


y-1  v-z 


T > r 
1 2 


_ Z_  cos  p d. 
r^=0  ^ ^^'*■^1 


y-1  2 ' 


1 (,l)l=+y-2^  . . . 


, , y-1  ijiw. 

y-1  n.  e 


sio 


f=0 


ij5. 


k-i-s 
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iU-i)v 


V-i  ] 

-rio  4+ri  ^ 


if  -pj^  is  a root  and  p is  odd. 


= (-l)^'*’^"^x^  . . . j,io  ^ ^k+r"  otherwise, 


.(P) 


'~  = "‘^k  ^ ^ 

& 

o 

= (-1)^^^  otherwise. 

From  these  relations,  it  is  clear  that  a^^|  and  "wiH  show 

regular  geometric  progression  behavior  for  large  consecutive  values 
of  p,  while  all  the  coefficients  between  will  be  irregular  because 
of  the  trigonometric  factors  involved. 

For  small  multiplicities,  }J  , the  trigonometric  relations  could 
be  used  to  determine  the  by  finding  cos  and  cos  (P+l)^jj.^j, 

and,  hence,  cos  For  larger  multiplicities,  the  resultant  procedure 

mentioned  in  Chapter  I might  be  employed;  p^,  of  course,  would  be  found 
from 


(3.26) 


Pk  - 


(P+1)  (p) 

\+U-l  \-i 

“(p7 

®k-l  \+i/-i 


CHAPTEE  IV 


SUIvjMARI  AMD  COMCLUSIONS 

Classical  root-squaxing  and  more  recent  root-cubing  axe  found  to 
be  special  cases  of  more  general  root-powering  as  expressed  in  the 
fundamental  theorem.  The  coefficients  of  the  powered  equations  are 
derived  foimally  in  three  different  ways;  as  convolution  products, 

'by  means  of  circulants  and  in  terms  of  Newton's  symmetric  functions. 

A fourth  procedure,  powering  of  the  companion  matrix,  is  also  discussed. 

Although  the  convolution  representation  of  the  coefficients  is 
formally  concise,  its  implementation  for  high  powers  in  numerical 
examples  becomes  quite  tedious.  Also,  it  lacks  the  desirable  property 
of  yielding  consecutive  powers  in  any  iterative  sense.  That  is,  each 
new  power  requires  a completely  separate  and  independent  treatment. 

For  practical  utilization  of  the  circulant  method  for  high-degree 
equations,  computer  programming  of  both  the  recurrence  relations  and 
the  initial  determinant  calculations  would  seem  to  be  necessary.  How- 
ever, once  started,  the  continuing  calculation  of  coefficients  for  suc- 
cessively higher  powers  is  simple  and  iterative. 

The  Newton  symmetric  function  procedure  offers  no  special  dif- 
ficulties for  a desk  calculator,  except,  possibly,  in  tabulation  of 
S-values  for  later  use.  It  appears  to  this  writer  that  this  method 
might  well  be  further  exploited. 

I>iatrix-powering  offers  particularly  efficient  and  easy  methods  of 
powering,  but  entails  the  problem  of  determining  the  powered  character- 


istic  equations.  Further  investigations  of  these  powered  matrices  with 
respect  to  their  characteristic  equations  or  their  characteristic  roons 
might  lead  to  useful  results.  The  traces  of  the  successively  powered 
matrices  offer  another  way  of  finding  the  Newton  symmetric  functions 
needed  for  that  method.  That  is, 

(4.1)  tr(A)  = , tr(A^)  = , . . . , tr(A^)  = . 

These  S's  can  then  be  used  in  (2.27)  or  (2.28)  to  find  the  coefficients 
of  the  powered  equations. 

Essentially,  from  the  consecutively  powered  equations,  the  actual 
roots  are  estimated.  Besides  the  obvious  advantage  over  root-squaring 
of  being  able  to  determine  the  unpowered  complex  roots  by  elementary 
methods,  there  is  the  possibility  of  finding  successive  approximations 
to  the  roots.  That  is,  the  roots  can  be  estimated  at  each  successive 
powering  and  their  convergence  used  as  a criterion  for  their  precision 
and  for  the  termination  of  the  powering  process. 


APPENDIX 


NUMERICAL  EXAMPLES 


The  roots  of  the  equation, 

(A.l)  (2x+l)(x+2)(x-5)  = 2x^  - 5x^  - 23x  - 10  = 0 


are  powered  by  the  usual  root-squaring  method  and  by  the 

four  methods 

developed  in  Chapter  II. 

Using  these  results,  the  roots 

are  approxi- 

mated  by  the  methods  of  Chapter  III. 

Root-squaring  method 

Power 

1 

2 

-5 

-25 

-10 

2 

If 

117 

429 

100 

(A.2)  If 

16 

10257 

160641 

1(4) 

8 

256 

10007(4) 

25600(6) 

1(8) 

16 

65536 

10001(12) 

65534(16) 

1(16) 

The  notation,  (n). 

indicates  the  power-of-ten  factor.  The  equation 

whose  roots  are  -(■ 

1 \l6 

> 

-(-2)^^  and  -(5)^^ 

is,  then. 

(A. 5)  65556x^  + 10001(12 )x^  + 6553^ (l6)x  + l(l6)  = 0 , 

or  the  equation  whose  roots  are  (-2)^^  and  (5)^^  is 

(A.4)  65556x^  - 10001(12)x^  + 6553^^(l6)x  - l(l6)  = 0. 

Convolution 

We  find  the  coefficients  of  the  equation  whose  roots  are  the 
fourth  powers  of  the  roots  of  (A.l)  using  the  relations 


^3 


2 3 

oc^  = i,  oc^  = -1,  = -i  and  performing  the  convolution  products  in 


the  order, 


}‘W 


(A.5) 

K h ^2  s \ s ^ ^ 

2 -5i  23  lOi 

2 3 -23  10 

4 lO(l-i)  25i  135 (1+i)  -529  230 (l-i)  lOOi 

2 3i 23  -lOi  

8 20  142  625  -1833  3315  -9667  4790  -2300  1000 

2 -3 -23  -10 

16  0 0 0 -10257  0 0 0 i6o64i  0 0 0 -i(4) 


The  desired  equation  is 


(A.  6)  l6x^  - I0257x^  + l6o64lx  - 10000  = 0 . 


Circulant 

The  six  hy  six  circulant,  from  which  all  powered  equations  for 
this  case  are  derived,  is 


-5x 

-23x 

-10 

0 

o' 

2x 

-5x 

-23 

-10 

0 

0 

2x 

-5 

-23 

-10 

0 

0 

2 

-5 

-23 

-10 

0 

0 

2 

-5 

-23 

-10 

0 

0 

si 

The  half-sequences  and  their  iteration  formulae  are  (A. 8) 


P = r+1  p = r 

; It,  5,  6 j = 5,  6 [ = 2 [j  4,  5,  6 1 

j;  3,  5,  6|  = 10  [;  1,  5,  6 
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1; 

5, 

61 

= 

10  1; 

Ij 

61 

li 

3, 

51 

= 

10  1; 

1, 

51 

1 s 

2, 

5, 

6l 

= 

2 1; 

3, 

5, 

6^ 

+ 

1 i 

2, 

61 

= 

-2  1; 

2, 

5, 

61 

\i 

2, 

51 

= 

“2  1; 

2, 

6i 

1 } 

1, 

5, 

61 

= 

2 1; 

2, 

5, 

61 

+ 

\i 

1, 

61 

= 

“2  1; 

1, 

5, 

61 

\i 

1, 

5i 

“2  1; 

1, 

6J 

‘ J 

2, 

5, 

61 

sz 

-10  1; 

1, 

3, 

6l 

1 J 

2, 

3, 

51 

= 

-10  i; 

1, 

3, 

5i 

1 i 

2, 

3, 

4f 

= 

-10  I; 

1, 

3, 

4l 

1 i 

1, 

3, 

61 

zz 

-10  1; 

1, 

2, 

6l 

* } 

1, 

3, 

5^ 

ss 

-10  I; 

1, 

2, 

5| 

1 > 

1, 

3, 

4r 

= 

-10  |; 

1, 

2, 

4( 

1, 

2, 

6r 

= 

-23  l; 

1> 

2, 

6r 

+ 

1 J 

1, 

2, 

5r 

= 

-5  l; 

1. 

2, 

6t 

+ 

1 J 

1, 

2, 

^Ir- 

= 

2 1; 

1, 

2, 

61 

1; 

1, 

2, 

31 

s 

10  1; 

1, 

2, 

31 

Listing,  in  the  same  order  as  (A. 8), the  values  of  these  with  their 


cofactors 

gives : 

Cofactor 

p=6 

P=7 

CO 

II 

p< 

S 

8 

16 

32 

2 

20Xp 

-710 

-6250 

-65910 

-92x^ 

100 

1420 

12500 

40x  p 

-40 

-200 

-2840 

142x; 

-1733 

-15795 

-164093 

-2?ox; 

270 

5466 

31590 

lOOXp 

-92 

-540 

-6952 

625x^^ 

-625 

-6591 

-64545 

-1733X 

142 

1250 

13182 

710x 

-20 

-284 

-2500 

958x 

-2300 

47900 

-966700 

-46Qx 

-500 

13500 

-285500 

2C0x 

200 

-4600 

95800 

2855X 

-4790 

96670 

-1957910 

-I55OX 

-1350 

28550 

-579150 

5OQX 

460 

-9580 

193340 

-9667X 

-9667 

193791 

-3878045 

479OX 

-2855 

57915 

-1162295 

-250OX 

958 

-19334 

387582 

1000 

-1000 

-l(^) 

-1(5) 

. . p=l4 

p=15 

P=l6 

2048 

4096 

8192 

-64939(6) 

12985(6) 

-25984(5) 

-16234(7) 

32464(6) 

-64952(5) 

-64954(6) 

12987(6) 

-25970(5) 

-64934(7) 

■12987(7) 

-25970(6) 

-16254(8) 

32469(6) 

-64928(6) 

-64936(7) 

129867(6) 

-25974(6) 

-64936(8) 

12986(7) 

-25974(7) 

-16234(9) 

52467(8) 

-65938(7) 

-64936(8) 

12987(8) 

-25973(7) 

-62060(9) 

-18618(9) 

62060(8) 

-12412(10) 

-37236(9) 

12412(9) 

-24824(10) 

-74472(9) 

24824(9) 

12412(11) 

37236(10) 

-12412(10) 

24824(11) 

74472(10) 

-24824(10) 

49648(11) 

14894(11) 
-49648(10)  1 

-24824(12) 

-74472(11) 

24824(11) 

-49648(12) 

-14894(12) 

49648(11) 

-99296(12) 

-29789(12) 

99296(11) 

-1(11) 

-1(12) 

-1(13) 
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The  sum  of  the  binary  products  of  a term  in  the  cofactor  column  times 

the  corresponding  term  under  an  indicated  p gives  the  equation  whose 
'til 

roots  are  the  p^  powers  of  (A.l).  For  p = 15,  and  l6  these  axe, 

respectively; 

16384x^  - 99999 (9 )x^  + l6584(l4)x  - l(lif)  = 0 
(A.9)  '^2160^  - 99999(10 )x^  , 32768(l5)x  - 1(15)  = 0 

65556x^  - 10001(12)x^  + 65554(l6)x  - l(l6)  = 0 , 


Matrix -powering 

We  first  divide  (A.l)  through  by  2 to  obtain 
(A.IO)  x^  - 2.5x^  - II.5X  -5=0. 

For  ease  in  tabulating,  the  order  of  the  rows  of  the  companion  matrix 
of  (A. 10)  is  reversed  and  the  new  rows  for  each  higher  power  added  at 
the  bottom. 

Power 


0 

1 

0 

1 

0 

0 

1 

2.5 

11.5 

5 

2 

17.75 

35.75 

12.5 

5 

78.125 

216.625 

88.75 

4 

4il.9575 

987.1875 

390.625 

For  the  equation  whose  roots  are  the 

powers 

we  transform  the  matrix 

(A.ll) 

k 

(411.9375 

987.1875 

PT  = 

78.125 

216.625 

17.75 

35.75 

of  the  roots  of  (A.lO), 


by  the  method  of  Danilevsky  [8,  pp.  166-I76]. 


390. 625 \ 
88.75 
12.5 


kS 


411.9375 

987.1875 

390.625 

78.125 

216.625 

88.75 

17.75 

53.75 

12.5 

-.5259259 

-.3703704 

17-75 

- 107.25 

29.25 

25 

53.75 

- 55.8037 

6.418518 

8.518519 

12.5 

0 

1 

0 

-3112. 062~ 

748.3125 

731.25 

.2404555 

.”2349728 

-3112.062 

.03446268 

3.461149 

-.2008316 

7J+8.3125 

1 

0 

0 

731.25 

0 

1 

0 

641.063 

-10040.1 

625 

The  desired  equation  is. 

then. 

(A.12)  x^ 

- 641.063x^ 

+ 10040. lx  - 625  = 0 . 

This  agrees  to  five  significant  figures  to  the  exact  coefficients  from 
root-squaring  (A.2). 


Symmetric  functions 

/'  Here,  again,  we  use  the  equation  in  the  form  (A.IO)  and  find  the 
equation  whose  roots  are  the  powers  of  (A.lO).  For  this,  we  need 
the  values  of  up  to  = S^j  they  are  calculated  here  hy  (2.25). 

51  = -a^  = 2.5 

52  = -aj_S^  - 2ag  = 2.5  X 2.5  + 25  = 29.25 

S " "^1^2  “ ^2^1  " ^^5  " -^75 

S]^  = -a^^S^  - a^Sg  - a^Sj_  = 61^1.0625 

(A.  15)  S " “ ®’2^  “ ^^2  ""  5°92.969 

Sg  = 15689.02 
= 77996.99 

53  = 390881 
= 1952613 

S^Q  = 976661^9 


^7 


= 48826080 

S12  = 244144700 

The  coefficients,  A^,  of  the  powered  equation  axe  found  from 

these  hy  (2.28); 

A^  = = -641.0625 

(A.14)  A^  = 4(83  + A^S^)  = 10040.1 

A5  = -^(S^  + A^Sg  + AgSj^)  = -625.15  . 

Determination  of  roots 

+V> 

From  the  last  two  equations  of  (A. 9)  for  the  I5—  and  16-—  powers, 
we  calculate  the  roots  by  (5.6). 

(16)  (15) 

V A ^ - 10001(12)  X 52768  _ _ - 

^ ^(16)  (15)  - 65556  X 99999(10)  ~ 

o  ^ 

a(l6)  (15) 

X - ^ ^ = 65554(16)  X 99999(10)  _ 

^2  (16)  (15)  10001(12)  X 52768(15)  ~ " 1*9997 

®-2 

(16)  (15) 

..  ^ 5 ^2  1(16)  X 52768(15) 

“5  -(16)  ^(15)  " “ 65554(16)  X 1(15)  = " • 


Similarly,  using  the  first  two  equations  of  (A.  9)  for  the  14—  and  15^ 


powers j 


„ i 99^99(10)  X 16584  _ _ --- 

1 ~ 52768  X 99999(9)  " 

X = - 52768(15)  X 99999(9)  _ ^ ... 

2 99999(10)  X 16584(14)  " ” 

V i 1(15)  X 16584(14)  _ 

5 “ " 5276b(15)  X 1(14)  " " *5°°^  ^ 


hQ 


and  we  have  agreement  in  all  cases  to  at  least  three  significant  figures. 
For  an  example  with  repeated  and  complex  roots  we  use 

(A.15)  - ,6x^  - 13. 2x^  - 36x  - li-0  = 0 

whose  roots  are  known  to  be  -2,  ^(-6^81).  The  equations,  found  by 

the  use  of  Newton's  symmetric  functions,  for  the  9 — f 10^,  11—  and 
‘til 

12 — powers  of  the  roots  are,  respectively, 

x^  - 1953096x^  - 557T85(2)x^  - 283(8)x  - 262144(9)  = 0 

x^  - 9764625x^  - 9770955 (5 )x^  + 1(13 )x  + 1048576(10)  = 0 
(a.  16)  k ^ p 

x^  - 48829(3 )x^  + 42817(6)x  + 4/3  (14 )x  - 4194304(11)  = 0 

x^  - 2441458 (2 )x^  + 12632 (8 )x^  +.553(16)  + 16777216(12)  = 0 

Immediately,  the  second  coefficient  is  seen  to  form  a geometric  progres- 
sion, the  next  two  are  irregular  and  the  third  is  regular  as  always. 

Thus,  the  largest  root  is  given  by 


V = -48829(5) 
1 -9764625 


5.0006 


-2441458(2) 
1 “ -48829(3) 


5.00002 


as  two  successive  approximations.  The  remaining  three  roots  must  have 
equal  moduli,  this  given  by 


3 _ 16777216(12)  X 48829(5) 

^ " 24415B(2)  X 4194304(11)  ~ 7*99997 

or  p = 2.  Examination  of  the  signs  of  the  second  and  last  coefficients 
indicate  the  remaining  real  root  to  be  - V 7*99997  = -2.  The  ampli- 
tudes of  the  remaining  two  conjugate  imaginary  roots  are  estimated  by 
use  of  trigonometric  relations. 


(A.17) 
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(m) 

= “p“  [(-ir  + 2 cos  m ^ ] 


a. 


\in 


for  m = p-1  = 10,  p = 11,  p+1  =12.  ^ is  not  used  in  this  case, 

"because  of  its  loss  of  significant  figures  due  to  round-off  error. 
With  some  elementary  manipulation,  we  obtain  * 


a(^) 

^2  ^ ^2 


(A.18)  cos  ^ = 


H. 


m * ~7iw  * ^ 


12 


pa- 


(11) 


12 


Using  the  data  from  (A.16)  and  p = 2 in  (A.18)  gives  cos  ^ = -.6OOO7. 


2 

P 

cos  <p 


Estimated  value 

5.00002 
-2.00000 
2.00000 
- .60007 


Actual  value 

5 

-2 

2 

-.6 
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